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ABSTRACT 



Stellar magnetic activity is a source of noise in the study of the transits of extrasolar planets. It induces flux variations which affect sig- 
nificantly the transit depth determination and the derivations of planetary and stellar parameters. Furthermore, the colour dependence 
of stellar activity may significantly influence the characterization of planetary atmospheres. Here we present a systematic approach 
to quantify the corresponding stellar flux variations as a function of wavelength bands. We consider a star with spots covering a 
given fraction of its disc and model the variability in the UBVRIJHK photometric system and in the SpitzerflRAC wavebands for 
dwarf stars from G to M spectral types. We compare activity-induced flux variations in different passbands with planetary transits 
and quantify how they affect the determination of the planetary radius and the analysis of the transmission spectroscopy in the study 
of planetary atmospheres. We suggest that the monitoring of the systems by using broad band photometry, from visible to infrared, 
helps to constraining activity effects. The ratio of the relative variations of the stellar fluxes in short wavelength optical bands (e.g., 
U or B) to near infrared ones (e.g., J or K) can be used to distinguish starspot brightness dips from planetary transits in a stellar light 
curve. In addition to the perturbations in the measurement of the planetary radius, we find that starspots can affect the determination 
of the relative semimajor axis and the inclination of the planetary orbit which have a significant impact on the derivation of the stellar 
density from the transit light curves. 

Key words. Infrared: stars - Techniques: photometric - Stars: activity - Stars: solar-type - (Stars:) starspots - (Stars): planetary systems 



1. Introduction 

The microvariability of active stars can hamper the measure of 
planetary transits and this can be a problem for targets with a 
solar-type or slightly higher magnetic activity level, i.e. main- 
sequence stars from late F to late G-type, or even worse for 
dM stars. The transit of a large group of spots on the solar sur- 
face causes a relative de crease of the optical flux of ~3.3 xlO 3 
dFrohlich & Leanl [2004) and even in the absence of sunspots, 
the solar flux exhibits a modulation due to bright structures, 
such as photospheric faculae. For comparison, the transit of a 
planet across the solar disc would produce a relative flux varia- 
tion of the order of 10~ 4 and 10~ 2 for Earth and Jupiter, respec- 
tively. Therefore, the search for planetary transits in the presence 
of stellar activity requires techniques to filter out or fit the mi- 
crovariability of the star. They are based on the observations of 
the flux variations outside transits along time intervals of at least 
two or three s t ellar rotations (cf., e.g., Ai grain & Irwinl 2004 
Moutou et al., 2005] iBonomo & Lanzal l2008t iBonomo etal 



2009). 



The detection of a transiting exoplanet is favoured at infrared 
wavelengths, where activity variations are smaller, especially 
when surveying dM-type stars. Moreover, their smaller size pro- 
vides a more advantageous ratio between the planet and the host 
star radii, implying a greater transit depth for a given planetary 
radius, and sugge sts that for dM stars it could be easy to detect 
low mass plane ts (Butle ret all l2004t INutzman & Charbonneaul 
2008; Deming et al.Ll2009h . However a large fraction of dM stars 



shows a significant level of magnetic activity which may be a 
crucial problem for the detection of transiting planets. 

In the infrared bands, the ratio of the planet to the star emis- 
sion increases remarkably, thus we can directly detect the ra- 
diation coming from the planet, allowing us to characterize the 
atmosphere during the occultation phase, i.e. through the analy- 
sis of the reflected planetary s pectrum (emission sp ectroscopy, 
see ICharbonneau et al.l 1 1 9991: ICam eron et al. 1999) or during 
the planetary transit in front of the star, i.e. by analyzing the 
transmission spectrum of the planetary atmosphere (transmis- 
sion spectroscopy, see ISeager & Sasselovl IToOOt [grown] [2001; 
iHubbard et al.ll2001t ICharbonneau et aLll2002l) . The analysis of 
infrared spectra is crucial for transmission spectroscopy, since 
the infrared domain includes the most important spectral sig- 
natures of the main atmosp heric molecu l ae even if their de- 
tection is still controversial ( Tine tti et al.L l2007at ISwain et all 
1 20081: iDesert et all [2009t ISing et all 120091: lAgol et all 12 010; 
Desert et all 1201 laHPalle et all 1201 lllGibson et al.Ll201lh ! 

The observation of stars hosting planets in the optical 
and near-infrared wavebands is the main scientific objec- 
tive of a number of grou nd-based projects (e.g. MEarth, see 
INutzman & Charbonneaul [2008b . It is par t of the program of 
the JWST (James Webb Space Telescope. iGardner et al.l 12006) 
and will be the focus of the proposed ES A mission, ECh Q}\ 
(Exoplanet Characterization Observatory, iTessenvil l2010h . a 



The EChO mission has been selected for the assessment phase for 
the medium-class mission, to be launched in the period 2020 - 2022 by 
ESA. 
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space-borne telescope that will observe the visible-infrared spec- 
tra (0.4 - 16 yum) to characterize the physical properties of ex- 
oplanet atmospheres, searching for molecular spectral features 
and bio-markers. A significant part of the mission will be de- 
voted to the observation of the atmospheres of Super-Earth plan- 
ets orbiting M dwarf stars for the aforementioned reasons. 

The most widely used approach to better constrain the plan- 
etary atmospheres is that of measuring the transit depth in dif- 
ferent passbands observed simultaneously. The selective absorp- 
tion of the stellar radiation by the molecular species, present 
in the planetary atmosphere, may produce the variation of the 
apparent radius of the planet versus the central wavelength of 
the pas sb and and can be us ed to detect them (e.g.. lDesert et all 
12011 alibi ICroll et all 1201 lb . The presence of stellar magnetic 
variability, due to cool spots, bright faculae, or magnetic sur- 
face inhomogeneities in general, can modify the transit depth 
and this can have a significant impact on the derivatio n of the de 



pende nce of the planetary radius on the wav elength (Pont et al 



2008: ICzesla et all l2009t ISing et all 12009b lAgol et all 12010b 
Bert a et alll201 UlSing et all hoi lab; Des ert et alll201 lal) . The 



issue is to understand if the apparent radius variations are due to 
molecular species in the planetary atmospheres or to stellar ac- 
tivity whose effect is also expected to depend on the bandpass. 
In this respect, the quantification of activity-induced variations, 
as functions of wavelength, may be crucial for the transmission 
spectroscopy and the characterization of the planet atmosphere 
and composition. These effects are crucial above all for active 
sources, such as HD 189733, but the issue rises even for less 
active sources, such as GJ 1214. 

Since the disc overy of its planetary companion 
dBouchv et all 120051) . HD 189733 has been th e subject of 

op 

2006t IWinn etaLl 120071: iPont et all 120071 12008b ISing et al 
2011b) an d in the infrared waveban ds dBeaulieu et al 



intense observations bot h in the op t ical bands dBakos et all 

~ H I2008P 
eaulieu e 
_s a K -tyr, 

a strong chromospheric activity dwright et all |2C 



2008; 



Sing et all 120091: iDesert et all 12011 ah. It is a K - type s tar with 

' ' 2004 which 

Toduces a quasi- periodic optical flux variation of ~1.3% 



pr 

dWinn et ail [2007) due to the rotation of a spotted stellar 
surface. The first accurate deter mination of the p lanetary 
system parameters is carried out by Bak os et al. (2006) through 
BVRI photometry; neglecting the effects of stellar variability, 
their planet-t o-star radius ratio ( 0.156 + 0.004) resu lts smaller 
than t hat of iBouchv et all 120051 (0.172 + 0.003fl IWinn et all 
(200 7]) obtain radius measurements compatible with those 
of iBakos et all (120061). agai n neg lecting the effects of stellar 
variability. E ven IPont et al.l d2007l) obtain results compatible 
with those of IBakos et a l. (2006), in this case correcting their 
HST/ 'ACS observations for the stellar variability; moreover these 
authors expect that the starspots effect be reduced in the infrared 
because of a lower spot contrast. Through //5T/NIC MOS transit 
observations, correcting for the unocculted starspots, Sing et al. 
d2009t) find a planet-to-star radius ratio of 0.15464 + 0.00051 
and 0.15496 + 0.00028 at 1.66 and 1.87 yum, res pectively, in 
agreem ent with the planet-to-star radius ratio of IBakos et all 
(2006j). In a multiwave l ength s et of HST/STIS optical transit 
light curves, ISing et all (1201 lbl) observe the typical signature 
of an occulted spot and, considering the effects of unocculted 



2 Throughout this section, we will present the results achieved in the 
literature by comparing the planet-to-star radius ratio of each source, so 
as to avoid any inconsistency due to the use of different stellar radii and 
to match quantities directly related to the observations. Where needed, 
we will compute the R p /R-„ ratio by using the stellar radius adopted in 
the relevant paper. 



spots, estimate a spot correction to the transmitted spectrum 
of about 0.00292 + 0.00113 as a fraction of the transit depth 
in the 320 - 375 nm passband, greater than in the 575 - 625 
nm passband where it is negligible. Moreover, they clearly 
demonstrate the dependence of the planet-to-star radius ratio 
on the out-of-transit stellar flux showing that the transit is 
deeper when the star is fainter as expected as a consequence 
of unocculted spots present on the disc of the star during 
the transit (se e Sect . 12.2 1. In Spitzer/IRAC wavebands, both 
IBeaulieu etail (l2008h and besert et al.l (1201 lah detect variations 
in the planet-to-star radius ratio. The former authors find that 
th e infrared transi t depth is smaller than in the optical (cf. Fig. 4 
of IBeaulieu et all [2008) and the effect of stellar spots is to 
increase this depth. By correcting for this effect, they find that 
the transit depth is shallower by about 0.19%, at 3.6 yum, and 
0.18%, at 5.8 yum. In the 5/?/fzer/IRAC's 3.6 yum band, the latter 
authors observe a greater Rp/R* ratio in low brightness periods 
because of the presence of starspots, implying that the apparent 
planet-to-star radius ratio varies with stellar brightness from 



0.15566!™^' to 0.1545 ± 0.0003 dDesert et all 12009) from 



low to high brightness periods, respectively. 

A less active M-type star which shows a rotational flux 
modulation of ~2% with a period of ~50 days is GJ 1214 
dCharbonneau et all 120091) . Considering the error bars, the 
planet-to star radius ratio measurements through optical and 
near-infrare d observ ations (ISada et all 120101: iBerta et ail 1201 ll 
ICarter et all [2011) are all consistent with the results of 



Charbonn eau et al.ll2009l (0. 1162 + 00067 ), obtained in the op- 
tical bands. Furthermore IBerta et al.l d201 lb detect variations in 
the planet-to-star radius ratio at different epochs and assess the 
expected variation of this ratio due to stellar variability in opti- 
cal wavelengths (see their Fig. 8). Th ey state that this e ffect may 
be neglected in Spitzer/IRAC bands. ICroll et aD (1201 ll) observe 
a deeper transit in the band than in the J band and, after tak- 
ing into account to effects of the starspots, explain this effect as 
due to a source of absorp tion in the 2.2 - 2.4 pm range, possibly 
associated with methane (Mi ller-Ricci & Fortnevl [20101) . 

A remarkable study of the starsp o t effec ts on the exoplanet 
sizes is c arried out by |Czesla et all (l2009h on the active star 
CoRoT-2 (lLanza et all 120091: iHuber et alll2010l) . They point out 
the importance of the normalization of the transit profile to a 
common reference level to make different transits comparable 
with each other and show that the relative light loss during a 
transit is significantly correlated with the out-of-transit flux. The 
correlation is compatible with a distribution of starspots occulted 
during the transit that has a larger covering factor when the star is 
more active, i.e., when the out-of-transit flux is lower. A method 
to derive the unperturbed transit profile is proposed by extrapo- 
lating the observed profiles towards their lower flux limit. A best 
fit to this profile gives a planet radius about 3% larger than that 
of lAlonso et al. (2008), which is estimated without considering 
stellar activity in the CoRoT white bandpass (350 - 1 100 nm). 

The above examples point out how stellar magnetic activity 
has a significant impact on the transit depth and in the conse- 
quent derivation of the planet radius and density, both in the op- 
tical and in the infrared wavebands, thus leading in several cases 
to contrasting results. As shown, the study of stellar variabil- 
ity is limited to its detection during individual real observations, 
without a systematic theoretical treatment to better characterize 
its impact through the analysis of stellar colours. This deficiency 
may lead to create "ad hoc" functiona l relations for each s pecific 
situation. (For instance, the Eq. (4) in lCarter et al.l (1201 ll) repre- 
sents the slight brightness increasing due to the starspot occul- 
tation during a planetary transit and holds only when a circular 
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Fig. 1. The transmittance of the filter set adopted in this work. 
Top: the Johnson-Cousins-Glass photometric system. Bottom: 
the Spitzer/IRAC passbands. 



spot, with the same planetary size, is considered.) A theoreti- 
cal characterization of transit light curves may help the correct 
definition of planetary and stellar parameters. 

In this study, we present a simple approach to better con- 
strain the spot contribution in the optical and infrared transit 
light curves, thus quantify the flux variations induced by stellar 
magnetic activity as a function of the wavelength passband, and 
investigate the starspot effects on the colours of stars. Moreover 
we present a unified theoretical approach in order to better con- 
strain the impact of both occulted and unocculted starspots on 
the determination of planetary and stellar parameters. Therefore 
in the context of the transmission spectroscopy, we examine the 
dependence of the apparent planetary radius on the magnetic ac- 
tivity level and wavelength. 



2. Method 

Our primary aim is to synthesize the flux variation in solar-type 
stars when active regions (dark spots in the considered case) 
stand on the stellar surface, to understand the effect induced on 
the light curve of a planetary transit. We focus on the optical 
and infrared bands, in the Johnson-Cousins-Glass photometric 
system (UBVRJJHK bands) and in the Spitzer/IRAC wavebands 
centered at 3.6 pm, 4.0 pm, 5.8 pm and 8.0 pm. 

We use the sta ndard BaSeL library of syn theti c stellar spec- 
tra, wo rked out bv iLei eune et al.l (1 1 997l 1 1 998h and lWestera et al.1 
d2002l) . Th is library is based on the gr id s of model atmosp heres 
spectra of iBessell et al.l d 19891 Il99 lh . iFluks et all d 19941) and 
Kuruczl d 19951) and covers a wide range of parameters: effective 
temperature 2000 K < T eS < 50000 K, gravity -1.02 < log g < 
+5.5, metallicity -5.0 < [M/H] < +1.0, wavelength 9.1 nm < A 
< 160000 nm. The stellar synthetic spectra have a wavelength 
resolution of 1 nm in the ultraviolet, 2 nm in the visible, 5-10 
nm in the near infrared, 20 - 40 nm in the Spitzer/IRAC bands. 



To derive the stellar flux in the Johnson-Cousins-Glass pho- 
tometric system and in the Spitzer/IRAC bands, we convolve the 
stellar flux distribution derived from the BaSeL libraries with the 
trasmission curves of the filters shown in Fig. [T| We choose as 
UBVRI filter passbands those ad opted b y Bessell dl990[) . as JHK 
filter passbands those adopted by Bessell & Brett ( 1988), and the 
Spitzer/IRAC filter passbands as reported in the Spitzer Science 
Center web sit^|. Both synthetic spectra and filter trasmission 
curves are interpolated in wavelength at the desired resolution. 

2. 1 . Modelling stellar activity 

We assume a temperature difference between the star and the 
spot, AT = T- k — T s , where and T s represent the effective tem- 
peratures of the unperturbed photosphere of the star (e.g. with 
no stellar activity ) and of the spot, re spectively. Following the 
same approach of iMarino et al . (1999), the measured flux of the 
system "star + spots", F* +s , corresponding to the specified tem- 
perature difference, is 



(1) 



where F+ and F s are the unperturbed and the spotted photo- 
spheric fluxes, respectively, / is the spot filling factor, or cover- 
ing factor, which represents the stellar disc fraction covered by 
spots, and A the mean wavelength of the passband. The relative 
variation of the stellar flux produced by starspots with respect to 
the unperturbed photosphere is: 



AF 
~~F 



F^X) - F^(A) 
FAX) 



= A A f, 



(2) 



where the contrast A ) = 1 -FJA)/Fj,(A) is a function of t he pass - 
band. Simil ar approache s may be found inlDesert et al.l (1201 lal) : 
iBerta et al.1 d2011l) and lCarter etail (1201 ll) . 

In the computation of the stellar flux, we assume that the stel- 
lar and the spot intensities vary according to the same quadratic 
limb-darkening law, viz.: 



I(A,p) = I(A, 1)[1 - mu(1 -Iti- «z»(l - AO 2 ]. 



(3) 



where I(A,p) is the specific intensity at wavelength A, p = cos 6, 
with 6 the angle between the normal to the stellar surface and the 
line of sight, and «y and u^a are the limb-darkening coefficients 
(LDCs) at the given wavelength. Given the difference in the ef- 
fective temperature, the LDCs of the unperturbed and the spotted 
photospheres are actually different, but in the modelling of the 
stellar and solar light curves they are assumed t o be the sam e to 
simplify the treatment (see, e.g., lLanza etafl I2003L 120041) . In 
our case, we have a maximum difference between the LDCs of 
~30% in the U passband when considering main-sequence stars 
with an effective temperature in the range 3700 - 6000 K and 
spots with a temperature lower by 1250 K. In view of the ad- 
vantages of assuming the same limb-darkening profile and the 
other uncertainties present in the problem, we shall neglect such 
a difference (cf. Sect. 12. 2\ . 

The effect of the wavelength depend ence of the limb da rken- 
ing is shown in the case of CoRoT-2 bv lCzesla et al.l l2009 using 
optical light curves at different wavelengths. The LDCs decrease 
from the U band to the infrared giving a more flat-bottomed tran- 
sit at longer wavelengths with shorter ingress and egress profiles 
(see discussion in Sect. 6 and F ig. 1 l. lSeager & Mallen-Ornelasl 
2003; iRichardson et all 120061) that can improve the determina- 
tion of the transit parameters. Thus even if we consider limb 



http://ssc.spitzer.caltech.edu/irac/calibrationfiles/spectralresponse/. 
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darkening in the derivation of the transit profiles, taking into ac- 
count that it decrea ses towards longer wavelengths (see Table 1 
in electronic form, IClaret et aUl 1995b . its impact will not be so 
crucial, sin ce our study is focu s ed primarily on infrared pass- 
bands (see iDeming et al.1 120071; ICarter et al.1 120081: iBerta et al 
l20ll . 

The stellar microvariability may be due also to the presence 
of faculae, bright photospheric regions surrounding spots. They 
exhibit centre-to-limb variation of their contrast and their bright- 
ness depends on size, position on the stellar disc and wavelength. 
In the case of the Sun, they can produce a relative in crease of the 
solar irradiance of ~10~ 3 ( Frohli ch & Leanl 120041) at the max- 
imum of the 11 -yr cycle. They are visible primarily in white 
light and the contras t to the unperturbed photosphere is reduced 
in the near-infrared (ISolanki & Unruhl Il998t iFrohlich & Leanl 
120041) . since it depends on A~ l dChapman & McGuirelll977h .ln 
the present study we neglected their contribution since we are 
interested in the infrared bands where the facular effect may be 
considered less important. Moreover, the facular contribution to 
the optical flux variations is negligible in th e case of stars sig- 
nificantly more active than the Sun (cf., e.g., lLanza et"al .1 120091 
and references therein). 

2.2. Effects of starspots on the apparent planetary radius 

Starspots affect in two ways the stellar flux variations in a tran- 
sit light curve. Spots that are not occulted by the planet produce 
a decrease of the out-of-transit flux of the star Foot that is the 
reference level to normalize the transit profile. Spots that are oc- 
culted during the transit produce a relative increase of the flux 
because the flux blocked by the planetary disc is lower than in 
the case of the unperturbed photosphere. The occultation of a 
starspot is therefore associated with a bump in the transit profile 
whose duration depends on the sp ot size and whose height de- 
pends on the spot contrast (e. g., [Pont et all 120071: [Sing et al., 
1201 lbl) . If we denote the flux during the transit as Fu, we have: 



FootW = F+(A) - AF S (A), 
and 

F u (A,n) = F 00J (A) - (I(A,fi)) nR 2 p , 



(4) 



(5) 



where AF S — A^/oF* is the flux decrease due to the unocculted 
spots having a filling factor fit (cf. Eq.|2]l, (/*(/!,//)) is the specific 
intensity averaged over the portion of the stellar disc occulted by 
the planet of radius R p . If we denote with /; the filling factor of 
the spots occulted during the transit when the planet's centre is 
at the disc position fi along the transit chord, then 

</(A/i)> = (l-WJA^V+MhiAf-t)) = (h(A,M))(l-fiA A ),(6) 

where (/*(/!,//)) is the average of the specific intensity of the 
unperturbed photosphere over the area occulted by the planet 
and (7 s (/l,ju)) is the same quantity for the spotted photosphere. 
Note that Eq. (O is valid because we have assumed that the 
spotted and unspotted photospheres have the same LDCs. If we 
define the relative transit profile at wavelength A as D(A, fi) = 
1 - F n (A,fi)/FooT(<l), after some simple algebra we find: 



(7) 



where R+ is the radius of the star and G(jj) is a function of the 
limb-darkening coefficients and the disc position // of the planet 
during the transit, but does not depend on the spot distribution 



and temperature. If we define the unperturbed transit profile as: 

(R \2 
jjfc 1 G(ji), the profile in the presence of spots is: 



D{A,ii)=D {A,ii) 



1 " Ma 



(8) 



If the ratio R^fR* £ 0.1 , we can assume that (I+(A,/j.)) 7*(/l,ju) 
(Mandel& Ago||200l) and the expression of the function G(ji) 
becomes particularly simple: 



1 - «u(l - fi) - u 2A (l -fif 

1 3 6 



(9) 



In Eq. (0 the filling factor of the occulted spots fi is in general a 
function of the position of the centre of the planetary disc along 
the transit chord, while the filling factor of the out-of-transit 
spots /o can be assumed to be constant because the visibility 
of those spots is modulated on timescales much longer than that 
of the transit, i.e., of the order of the stellar rotation period that is 
generally of several days. Moreover, when a geometric model is 
fitted to an observed transit profile, the bumps produced by oc- 
culted spots are excised from the data set to reduce their impact 
as much as possible. However, the capability of detecting the 
flux bump due to a particular spot depends on its size and con- 
trast as well as on the accuracy and cadence of the photometry. If 
there is a more or less continuous background of spots along the 
transit chord, they cannot be individually resolved and the tran- 
sit profile will appear to be symmetric and without any obvious 
indication of starspot perturbations. This is indeed the most dan- 
gerous case when fitting an observed transit profile. In this case, 
fi can be assumed to be constant during a given transit because 
otherwise we would have detected the individual spot bumps. If 
we assume that both /o and /; are much less than unity, we can 
develop the denominator of Eq. (0 in series and neglecting the 
second order terms in the filling factors, we find: 



D{A,fi) 



6V)|i-.u/i -./;.»|. 



(10) 



From Eq. ( TTOt . we see that the effects of the unocculted and oc- 
culted spots tend to compensate each other although the filling 
factor of the occulted spots may generally be higher because 
starspots tend to appear at low or intermediate latitudes in Sun- 
like stars. 

An alternative interpretation of Eqs. <[8J an d < fT0b is that spots 
change the apparent relative radius of the planet (Rp/Ri,)^ as de- 
rived from the fitting of the transit profile (see Sect. 12.3b as: 



Rr 



Rn 



l-fiA A 



Rn 



l-^/i-Zo) 



Ma \R*Jyl-foA A \R- 

Therefore, the apparent variation of the planetary radius is: 

AR P (A) 1 
-^--^(/i-/o), 



(11) 



(12) 



When radius measurements are available in two different pass- 
bands Ai and A2, we have: 



R r 



Rr 



1 



*7/ ^"-2 (A * - A ^-M-fo)- 



(13) 



If the spot temperature is known, this equation can be used to 
derive the effective filling factor / = /; - /o from the apparent 
radius measurements in two well-separated passbands, e.g., in 
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Table 1. Examined cases and their related real planetary systems. In the synthesized models the metallicity and gravity values are 
the solar ones. 





Reference 


Systems 




Real Systems 


Case 


Spectral 


T, 


T, 


AT 


Star 


Spectral 


Rp 


M„ 


Am 




Type 


(K) 


(K) 


(K) 


name 


Type 


(Rj) 


(My) 




1) 


Early-G 


5750 


4500 


1250 


CoRoT-12 


G2V 


1.44 ±0.13 


0.917 


0.018 


2) 


Early-K 


5250 


4000 


1250 


CoRoT-7 


KOV 


0.15 ±0.008 


0.0151 


3.1410" 4 


3) 


Mid-K 


4500 


3250 


1250 


HAT-P-20 


K7 


0.867 ± 0.033 


7.246 


0.016 


4) 


Eaiiy-M 


3750 


2500 


1250 


GJ436 


M2.5 


0.365 ± 0.02 


0.0737 


6.53-10- 3 



Table 2. Contrast coefficients for each reference system in the different passbands. 



Case 


Aa 


A B 


A v 


A R 


A, 


Aj 


A H 


A K 


A 3 .6 


A4.5 


A5.8 


A 8 .o 


1) 


0.938 


0.854 


0.759 


0.686 


0.616 


0.458 


0.316 


0.285 


0.248 


0.283 


0.262 


0.232 


2) 


0.958 


0.901 


0.833 


0.756 


0.643 


0.560 


0.456 


0.429 


0.392 


0.397 


0.373 


0.348 


3) 


0.960 


0.959 


0.946 


0.912 


0.779 


0.614 


0.605 


0.547 


0.495 


0.444 


0.496 


0.503 


4) 


0.994 


0.995 


0.992 


0.978 


0.924 


0.722 


0.733 


0.640 


0.542 


0.543 


0.656 


0.677 



the U and the K passbands. If the spot temperature is not known, 
it can be derived when a third measurement is available in a fur- 
ther passband, thereby removing the degeneracy between filling 
factors and contrast. However, the presence of a planetary atmo- 
sphere with a wavelength-dependent opacity can complicate the 
problem as we shall see in the final part of Sect. I3.2l in the case 
of HD 189733. 

2.3. Effects ofstarspots on the determination of the orbital 
parameters 

The perturbations of the transit profile due to starspots do not 
affect only the determination of the relative radius of the planet 
Rp/R-i, but also the other geometrical parameters of the system 
and the limb-darkening coefficients. Most of the pre vious studies 
have considered only the radius variation (cf., e. g., iBerta et ail 
1201 it ICarter et"ail 1201 but the i mpact on the other parame- 
ters can be significant. lCzesla et al.l (120091) explored the variation 
of the planetary radius and orbital inclination keeping fixed the 
semimajor axis of the orbit a, the stellar radius and the limb- 
darken ing coefficients u\ and 112 at the values of Alonso et al. 
d2008h . In this way, they found that the radius of the planet is in- 
creased by ~3% by fitting the lower envelope of the transit pro- 
files derived with their method. On the other hand, they found 
no significant correction in the inclination i of the orbital plane 

of the planet to the plane of the sky. 

A limitation of the approach bv lCzesla et al.l d2009l) is that by 
fixing the ratio a/R^, the duration of the transit fx fixes the incli- 
nation because tj - Vl - b 2 P OI b/(na/R+), where b - a/R^ cos i 
is the impact parameter and P or b the orbital period. Since the 
effects of the spots on the duration of the transit are generally 
small, the inclination is not changed by fitting their unperturbed 
transit profile. Therefore, to perform a proper estimate of the 
systematic effects produced by occulted and unocculted spots, 
we cannot fix any of the system parameters in fitting the transit 
profile. 

The unperturbed transit profile Do(A, fi) can be derived from 
the lower envelope of the observed transit profiles if we can 
guess the minimum spot filling factor of the unocculted spots 
fo m ; n and of the spots occulted during transits From the 

lower envelope of the observed transits D\ ow (A,/u), we get an es- 
timate of the unperturbed profile as: 



Since f\ ^ is a constant, the unperturbed profile does not show 
any bumps due to the occultation of individual spots and is sym- 
metric with respect to the mid of the transit. 

We fit the unperturbed pr ofile u sing an analytical model for 
the transit, e.g., that of iPall (120081) that has the advantage of 
providing analytic expressions for the derivatives of the tran- 
sit profile Do with respect to the transit parameters allowing us 
an efficient applica tion of the Leve nberg-Marquardt method to 
minimize the x 2 dPress et all Il992h . To minimize their mutual 
correlations, we choose as free parameters for our best fitting: 

= (2n/P olb )(a/R+)/ Vl - b 2 , b 2 , R p /R+, and the limb dark- 
ening parameters u+ _ u\ + M2 and u- = U\ — u-i (see IPall 
2008: [Alonso et aUl2008l for a justification of the independence 
of this parameter set). The orbital period and the mid transit 
epoch are kept fixed in the fitting because they can be derived 
from the observation of a sufficiently long sequence of transits. 
Specifically, the profile distortions due to resolved occulted spots 
may affect the timing of individual transits, but the spot pertur- 
bations can be assumed to be random so that they increase the 
statistical error of the timing. Therefore, a sufficiently long time 
series can be used to derive an accurate value of P 01 -b. We shall 
present an application of this method to estimate the impact of 
starspots on the fitting of the parameters of the system CoRoT-2 
in Sect. 133] 



D (A,fi) = D lov/ (A,fi) 



1 - A A f , 



(14) 



3. Results 

We focus our analysis on star temperature values ranging from 
the solar one (~5750 K) to that of an early M-type star and set the 
temp erature difference between the star and the spot to 1250 K 
(see|Berdvugina 2005). The synthetic spectra, both for the unper- 
turbed photosphere and the spots, were chosen for solar metal- 
licity and gravity, i.e. [M/H] G = 0.0 and log g G = 4.44. To apply 
our results to real cases, we selected four systems with plan- 
ets, i.e., CoRoT-12, wi th a transiting Jup iter-like planet orbiting 
around a G2V star (seelGillon et alll2010h : CoRoT-7, a KOV star 
dLeger et all l2009t , lOueloz et al.L 12009^ HAT-P-20. a K7V star 
(Bak os et all l2010h : and GJ 436, a M2.5 dwarf star ha rbouring 
a transiting Neptune-mass planet (Bu tler et all 120041) . We list 
their properties in Table Q] where we give, from the left to the 
right, the spectral type, the selected temperature for the quiet 
photosphere and spot for the reference stars, the difference be- 
tween photosphere and spot temperatures, the real system with 
the name of the star we will use for comparison, the stellar spec- 
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Fig. 2. Flux vs. Wavelength for all the cases in TableQ] Solid lines plot 
the spectral distributions of the stellar unperturbed photospheres and 
dotted lines the spectral distributions of the star + spot systems with/ = 
0.5. In the top left we present the first case, in the top right the second 
case, in the bottom left the third case and in the bottom right the fourth 
case. 

tral type, the planetary radius and mass, and the reduction in 
flux at mid planetary transit, expressed in magnitudes (obtained 
as the square of the ratio of the planet radius to the star radius), 
respectively. 

The contrast coefficients A,\, giving the spot-induced flux 
variations for each of the examined cases (see Sect. 12. Il l, are 
obtained by choosing the spot filling factor of the out-of-transit 
spots, fo, as rangin g from 0.01 to 0.5 of the stellar disc surface 
(Berdvugina, 2005) and are reported in Table|2] For G and early- 
K stars, the minimum flux perturbation occurs in the 8.0 yum 
SpitzerflRAC passband, while the maximum occurs in the U 
passband because the spot contrast increases monotonically to- 
wards shorter wavelengths. In the case of late-K and M stars, 
the minimum flux perturbation occurs in the 4.5 /urn and 3.6 fim 
SpitzerflRAC passband, respectively, due to molecular bands in 
the spectrum. These features are observed in the spectra of cooler 
stars, both in the optical and in the infrared, and this results in 
a non-monotonous spectral variation in correspondence of these 
bands. In this case, the unperturbed and spotted spectra do not 
behave just as power laws in the infrared domain as in the hotter 
cases. 

In Fig. [2] we plot the system flux vs. wavelength for all the 
cases, for a given filling factor of the out-of-transit spots. The 
solid curve represents the spectrum of the unperturbed photo- 
sphere and the dotted line the spectral distribution of the system 
"star + spots". We have plotted a case with an extreme filling 
factor (/o = 0.5) to make the effect of the spots clearly evident. 
The spot perturbation decreases with increasing wavelength, al- 
though the details of the variation depends on the stellar effective 
temperature. 

In Fig. [3] we plot the variations of the normalized ratios 
Fv/Fbaad vs. the filling factor fo, where Fband is computed in 
three specified passbands (I, K and 8.0 yum), with the purpose of 



Fig. 3. Normalized F v /F b:m( ] vs. Filling Factor f for all the cases in 
Table [T] in different passbands: I, K and 8.0 pm. Different curves plot 
different ratios: Fy/F^ Fy/F K and F v /F 80 , as labelled. In the top left 
we present the first case, in the top right the second case, in the bottom 
left the third case, and in the bottom right the fourth case. 

understanding to what extent the stellar magnetic activity may 
produce a variation in the colours of the star (this graph will be 
discussed in details in sect. 13. 11 1. 

These results suggest that the colour effect due to stellar ac- 
tivity has to be taken into account if we analyse the photomet- 
ric signals even of a moderately active star (comparable with or 
slightly more active than our Sun). The effects are both on stellar 
colours (discussed in Sect. 13. Il l and on the determination of the 
radius of the planet (discussed in Sect. I3.2I >. 

3.1. Influence of stellar activity on stellar colours and the 
characterization of a single light curve dip 

If we consider the first and the fourth cases in Table Q](i.e., the 
hottest and the coolest host stars), by fixing the filling factor fo 
to 1% (i.e., larger than in the solar case, but comparable to those 
observed in late-type moderately active stars), we can compute 
the relative variation of the system flux (AF/F) in the U and K- 
bands, to compare the influence of stellar activity in the optical 
and infrared passbands. 

For the first case (7* = 5750 K, AT = 1250 K) we derive the 
following variations in the U-band and in the K-band, respec- 
tively: 

AFu/Fu = 9 ■ IQ>- 3 and AF K /F K = 3 ■ 10~ 3 . (15) 

For the fourth case (T+ = 3750 K, AT = 1250 K) we de- 
rive the following values in the U-band and in the the K-band, 
respectively: 

AFu/Fu = 1 • 10~ 2 and AF K /F K = 6 • 10~ 3 . (16) 

The decrease in flux caused by a spot filling factor of 0.01 
is comparable with the decrease produced by the transit of a 
Jupiter-sized planet in front of the stellar disc. In fact we can 
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compare the obtained reductions in flux of the system due to 
stellar magnetic activity, using the information in Table [2] with 
the decrease in magnitude due to the transit of a planet, reported 
inTableffl 

Our results indicate that it is easier to identify the transit of 
a Jupiter-like planet by analysing the light curve of the system 
in the K-band, while it may be more problematic in the U-band. 
The detection of a smaller planet is more challenging, both in the 
U and in the K-band, since the photometric signal correspond- 
ing to the transit of such a planet is, as expected, smaller than the 
effect due to the magnetic activity of its parent star. Fortunately 
other properties of the variations, as timescales or colour depen- 
dence, may help in understanding the origin of the variations. 

The transit of a spot or a group of spots (or an active region in 
general) on the star surface induces a colour variation of the star. 
Since spots are cooler than the stellar photosphere, the colour 
indexes are increased. In principle, we can infer the activity of 
a star by monitoring the variation in the colour for some spec- 
ified wavebands. Star colours depend on the ratio between the 
fluxes in two selected passbands. For this purpose, we compute 
the ratio between the system fluxes in the optical and the in- 
frared passbands, for all the four template cases: Fy/Fj, Fy/F^ 
and Fy/F&.o. In Table [3] we report the variation of Fy/F^d in 
the different bands, by fixing the filling factor fo to 0. 1 and to 
0.01, respectively. The first column shows the filling factors and 
the enumeration of the cases, the following columns show the 
changes in Fy/Fi, in Fy/F^, and in Fy/F$x>, respectively. The 
behaviour of the star flux ratios vs. spot filling factor of the un- 
occulted spots are plotted in Fig. [3] where the ratios decrease as 
the spot filling factor increases. 

The transit of a planet without atmosphere induces a relative 
decrease of the stellar flux that depends on the mean wavelength 
of the passband owing to the wavelength dependence of the stel- 
lar limb darkening. For the present application, we consider the 
central transit (b — 0) of a small planet (/? p /7?+) ^0.1, with the 
quadratic limb-darkening law of Eq. (0. We derive the ratio % 
of the relative depths of the transit at the mean wavelengths of 
two different passbands, say U and K, 

^ ^ AFu/Fy = 1 ~ \ u ik ~ \uik 
AF K /F K i-L Uw -l U2u 

where u\u, U2u, u\k, and u-lk are the limb-darkening co- 
efficients in the U and K pass bands, deriv e d fro m, e. g., 
iDiaz-Cordoves et alJ (1 19951) and IClaret et al.1 (1 19951) . respec- 
tively. For the first case, we find % = 1.298, while for the fourth 
case, % = 1.197. This is significantly smaller than expected in 
the case of the variations induced by dark spots (cf. (Tl3T > and 
(TTrjl)). The flux in the K passband coming from the night side of 
the planet during the transit is neglected because the planet tem- 
perature is always lower than ~2000 K, even for the strongest 
irradiated objects, that is remarkably lower than the temperature 
of the starspots. 

Considering Eq. (0, we find that the ratio of the variations 
induced by starspots is: 



For a stellar effective temperature of 5750 K, this is always sig- 
nificantly greater than 1.3, except for very cool spots, i.e., Ar > 
2100 K. Indeed, sp ots so cool are gen erally not observed in ac- 
tive stars (see, e . g . . iBerdvu ginal l2005h . For the Sun as a star, Ar 
- 500 - 600 K because suns pot irradiance is d ominated by the 
penumbral regions (cf., e.g., lLanza etaTl 120041 and references 



Table 3. Variations in the colour of the system for all the cases. 
Top: computed for fo = 0. 1 . Bottom: computed for fo = 0.01 . 



/o = 0.1 


A(F V /F,) 


A(F V /F*) 


A(F V /F 8 . ) 


1) 


1.5% 


4.9% 


5.4% 


2) 


2.0% 


4.2% 


5.0% 


3) 


1.8% 


4.2% 


4.7% 


4) 


0.7% 


3.8% 


3.4% 


fo = 0.01 


A(F V /Fj) 


A(F V /F K ) 


A(F V /F 8 . ) 


1) 


0.14% 


0.48% 


0.53% 


2) 


0.19% 


0.41% 


0.49% 


3) 


0.17% 


0.40% 


0.45% 


4) 


0.07% 


0.35% 


0.32% 



therein). In the case of the star with = 3750 K, a spot tem- 
perature deficit Ar > 1500 K is required to have !R S comparable 
with that of a plane tary transit, which is also unlike ly for such a 
cool star (see, e.g.. lBerdvug"inal 120051: IZborill.120031) . 

In conclusion, even for an individual light curve dip, we are 
able to compare the variation of the depths in the optical and 
infrared passbands to discriminate between an effect induced by 
spot activity and a planetary transit, because the ratio % is in the 
range 1.7 - 3.0 in the former case (see the ratio of relations in 
(TT3T > and in (O, for a typical Ar = 1250 K), while it is around 
1.2 - 1.3 in the latter case. 

The timescales of the stellar flux variability and of the transit 
are generally quite different, since the form er may range from 
week s to months, as in the case of the Sun ( lLanza et al.L [2003. 
2007), while the latter has a duration ranging from a few hours 
to a few tens of hours. However, there are situations in which 
the timescales can be comparable, e.g., in the case of a young 
sun-like star with a rotation period of, say, 1-2 days transited 
by a planet at the distance of ~1 AU, whose transit duration is 
of the order of ~ 10 - 15 hours. In this case, a discrimination of 
the origin of the light dip based on timescales is unfeasible while 
our approach is still applicable if simultaneous multiwavelength 
observations are available. Note that our method can be applied 
in principle even to a single event. This can be very useful in the 
case of planets with orbital periods of several tens or hundreds 
of days. 

The photometric precision currently achievable from the 
ground in the near infrared is generally sufficient for an applica- 
tion of the present method, provided that the photome tric varia- 
tions are at least of 0.005-0.01 mag. lCroll et al.1 (120111) acquired 
photometry with a standard deviation of <x » 5 x 10~ 3 mag in the 
J (1 .25 i-im) and the K (2. 15 nm) bandpasses usin g the WIRCam 
at the 3.5-m CFHT, while Tofflemi re et al.l d201 ll) reached a pho- 
tometric stability of 3.9 x 10~ 3 mag in one hour and of 5.1 X 10~ 3 
mag in one night in the K s passband at 2.224 jjm observing sta- 
ble M-type stars. However, the situation is remarkably different 
in the mid and far infrared passbands. 

From space, the Warm Spitzer mission has reached a statisti- 
cal error of 450 ppm (parts per million) per m inute, observing GJ 
1214 for 12 sec at 3.6 /im (billon et all 1201 ll) . However, the IR 
observations are affected by systematic effects that may compro- 
mise the measurements, since they induce flux variations com- 
parable to the depth of a planetary transit. For example, the 3.6 
and 4.5 /.im Spitzer/IRA C channels (InSb detectors) are affected 
by the pixel-phase effect ([Reach et all 120051: ICharbonneauet all 
120051: iMorales-Calderon et all l2006t iKnutson et al l 120081) . due 
to the telescope jitter and intra-pixel variation of the sensitiv- 
ity of the detector, which may cause a flux peak-to-peak ampli- 
tude of up to -1% dBeaulieu et all l2008t iHebrard et all 120101) 
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or even higher (iDesert et all l2011bt iFressin etal ac- 
cording to the exposure time. The light curves of the 5.8 and 



tion of the radius and the orbital parameters will be illustrated 
in Sect. 13.31 The plots in Fig. [4] AR p /R p vs. wavelength, give us 



f =0.5 
A 



• CoRoT-12b 



f =0.5 



• CoRoT-7b 



f =0.5 



HAT-P-20b 



f =O.S 



■ GJ 436b 



10.00.1 

Wavelength (/-im) 



Fig. 4. Relative Variation of the apparent planetary radius vs. wave- 
length for all the examined planetary systems. In the top left we present 
the results for CoRoT-12b, in the top right for CoRoT-7b, in the bottom 
left for HAT-P-20b, and in the bottom right for GJ 436b, assuming var- 
ious levels of activity. The dashed curves, from the bottom to the top, 
correspond to increasing filling factor values, as indicated on the left of 
each plot. The solid curves refer to the activity levels as actually esti- 
mated for the specified systems. The filled points indicate the position 
of the observed planets, as described in the text. 

8.0 fim channels (Si:As de tectors) show non-linear trends vs . 
time, called the ramp effect (iDeming et alll2006tlKnutson et all 
2007a]), caus ed by the trappin g of electrons by the detector im- 
purities (see lAgol et al.ll2010l for further details), that produce 
mmag level flux variations in photometry. They can reach up to 
~10% for the lowest illuminated pixels over 3 3 hr (lAgol et all 
l201Ch . and increase over time (Machalek et al., 2008). The cor- 
rection of these instrumental systematics and artifacts is partic- 
ularly relevant for future infrared missions, such as JWST, since 
the same Si:As technology, present in the last two SpitzerflRAC 
cha nnels, will be adop ted for MIRI (Mid-Infrared Instrument, 
see iGardner et al.ll2006h . 

Anyhow, in this study we intend to emphasize the importance 
of a broad band monitoring of stars hosting planets from visible, 
at the shortest wavelength as possible, to IR. Simultaneous ob- 
servations, from optical to infrared wavebands, allow us to better 
understand the role of stellar magnetic activity in the analysis of 
a transiting planet light curve comparing the ratios % derived in 
eqs. ( ITTb and (1181 1. as stated before. 



3.2. Influence of stellar activity on planetary radius 
determination 



In this section, we apply the simple formulae of Sect. l2.2l to eval- 
uate the variation of the apparent planetary radius in the presence 
of stellar activity assuming that the other parameters are not af- 
fected. The effect of the activity on the simultaneous determina- 



Table 4. The relative variation in planetary radius computed at 
optical wavelength for CoRoT-12b, CoRoT-7b e GJ 436b, and 
in the R-band for HAT-P-20b, taking into account the adopted 
filling factor values reported in the text. AErr is the ratio of the 
computed AR p to the error bars quoted in the literature. 



Planet Name 


AR P /R„ 


AErr 


CoRoT-12b 


1.14 • 1(T 2 


13% 


CoRoT-7b 


1.67 • 1(T 2 


31% 


HAT-P-20b 


0.46 • 1(T 2 


12% 


GJ 436b 


0.50 • 10~ 2 


10% 



the overestimate of R p if the effect of unocculted starspots is ne- 
glected for our four reference cases. We find a relation between 
AR p /R p and the wavelength and between AR p /R p and the stellar 
activity, i.e. increasing the filling factor of the out-of-transit spots 
fo the dashed curves shift upwards and the variation decreases 
towards longer wavelengths. The filled circles on the plot indi- 
cate the correction on the observed planetary radius that we have 
to apply to get the true radius of the planet for the assumed spot 
temperature and filling factors (see below). 

In the top left panel of Fig.|4]we present the result for CoRoT- 
12. We assumed a filling factor of /o = 0.03 (the solid line in the 
panel) by inf erring the flu x variation from the CoRoT light curve 
of Fig. 1 in iGillon et all (1201 Oh and then derived the contrast 
coefficients from the first row of Table [2] 

The top right panel of Fig. [4] shows AR p /R p vs. wavelength 
for CoRoT-7. The spot filling factor value is fo = 0.04 by de- 
ducing it from the optical flux variability reported in lLanza et all 
d2010l) . 

For HAT-P-20, the outcomes are shown in the bottom left 
panel of Fig. |4] We set fo = 0.01 as derived by the magnitude 
variations o f the s tellar flux in the R-band shown in Fig. 1 in 
Bakos et al. (2010). The dependencies of AR p /R p on the wave- 
length and the stellar activity are clearly visible for this planet 
too, even though the curves tend to become flat. 

The bottom right panel of Fig. [4] shows the result for 
GJ 436. The star is a relatively quiet source showing photomet - 
ric variations of the order of the millimag (But ler et all [2004). 
iDemorv et all d2007l) find that the photometric irregularities in 
the light curve of GJ 436, detected with the EULER Telescope, 
are due to the presence of a unocculted starspot and they com- 
pute a spot filling factor of 1%, which is compatible with the 
dispersion in HARPS radial velocity. Thus we set fo = 0.01. 

In Table [4] we list the obtained results: from the left to the 
right the planet name, AR p /R p , and the relative variation in the 
radius uncertainty, AErr, which is the ratio of the computed AR p 
with respect of the uncertainties quoted in the literatur^3 (see 
also Table [TJ. By a simple check, we notice that the largest val- 
ues of AR p /R p and AErr occur for the most active star within our 
sample, i.e. CoRoT-7 with fo = 0.04, while the lowest values 
are found for the coolest and less active stars, i.e. HAT-P-20 and 
GJ 436 with fo = 0.01 in both cases. For CoRoT-7b, we com- 
pute an overestimate of the planetary radius, AR p , of about 180 
km at optical wavelengths (with a relative variation in the ra- 
dius uncertainty of about 31%), while in the K-band it is about 
92 km and even smaller in Spitzer/IRAC wavebands (at 3.6 /im, 
AR p =i 84 km and at 8.0 fim AR p 75 km). Such a decrease 



See site http://exoplanet.eu/index.php 
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Table 5. Corrected planetary radii in the optical and Spitzer wavebands for three notable examples: HD 209458b, HD 189733b and 
GJ436b. 



Planetary 


Bands 








Ref. 


Systems 






(R./) 


(Ry) 








fo = 0.004 


(fo = 0.03) 








V 


1.347 


(1.335) 


1.349±0.022 


(1) 


HD 2094 ^Hh 


3.6 yum 


1.351 


(1.347) 


1.351±0.006 


(2) 


4.5 yum 


1.355 


(1.350) 


1.355±0.008 


(2) 




5.8 yum 


1.386 


(1.382) 


1.387±0.007 


(2) 




8.0 yum 


1.381 


(1.377) 


1.381±0.005 


(2) 






Wo = 0.01) 


fo = 0.03 








V 


(1.200) 


1.189 


1.205±0.003 


(3) 


HD 189733b 


3.6 yum 


(1.182) 


1.177 


1.184±0.002 


(4) 


4.5 yum 


(1.191) 


1.186 


1.194±0.002 


(4) 




5.8 yum 


(1.184) 


1.179 


1.186+0.004 


(4) 




8.0 yum 


(1.182) 


1.177 


1.184±0.003 


(4) 






fo = 0.01 


(fo = 0.03) 






GJ 436b 


V 


0.400 


(0.396) 


402 +0 036 


(5) 




8.0 yum 


0.378 


(0.375) 


0.379±0.018 


(6) 



References. (1) iRnutson et alJ j2007bl) : (2) iBeaulieu et all d2010l) : (3) iPont et all ([2007); (4) iPesert et alJ d2009h : (5) iBean et alj d2008l) : (6) 
iDeming et all d2007h . 



in AR p as wavelength increases is exhibited by the other exam- 
ined planetary systems and could be mis-interpreted as evidence 
of a planetary atmosphere with a wavelength-dependent absorp- 
tion. Since AR p is proportional to the planetary radius, we ob- 
tain the largest AR p for CoRoT-12b (AR p ^ 1180 km at opti- 
cal wavelengths) with a Jupiter-like planet and the lowest for 
GJ 436b (ARp =s 130 km in the optical) with its Super-Earth 
planet. Nevertheless, all the radius increments calculated in this 
work fall inside the error bars quoted in the literature, thus no 
correction to the presently determined planetary parameters is 
required. We see that for each value of the filling factor, the 
curves seem to start nearly at the same level (i.e. AR p /R p 0.005 
for f = 0.01 and AR p /R p - 0.24 for f = 0.5), independently of 
the stellar photospheric temperature. Their trends keep on to be 
monotonic but their slopes decrease towards cooler photospheric 
temperatures (cf. GJ 436, in the bottom right panel of Fig. 0J. 
This indicates that the overestimate of the planetary radius de- 
pends on the wavelength, with the dependence being steeper for 
solar-type stars and flatter for later spectral types. 

Specifically considering the Spitzer passbands, we apply 
the method described above to two of the most studied plan- 
etary systems: HD 209458 and HD 189733. These systems 
have been the objects of intense studies in a wide range of 
wavelengths: f rom optical with the HST (iBrownet all [ 200 1 ; 
iKnutson et all l2007bl: IPont et all 120071: ISineetall 1201 lbl) to 
infrared with the Spit z er satellite (Richardson et al., 2006; 
Ehre nreich etafl 120071: iTinetfi et all l2007bt IBeaulieu et all 
2008. 120 10t IPesert et all 1201 lal) , also as targets for transmis- 



sion spectroscopy thanks to their brightness. 

In Table [5] we report the results and present the correspond- 
ing ones for GJ 436b for comparison. For all the cases, we have 
considered two different filling factor values of the unocculted 
spots: one is that argued or reported in the literature, while the 
other one is given for comparison (indicated in parenthesis in the 
Table). 

HD 209458 is an almost quiet solar-like GOV star and the 
calculations have been made considering two configurations: 
/o = 0.004 and fo = 0.03. The former fo value is estimated 
assuming a circular spot with an ave r age ra dius of 4.5 x 10 4 km 
as in the model proposed by ISilval d2003l) to fit the distortion 
of a single transit of HD 209458 observed by HST. HD 189733 



is an active early-K star with fo = 0.03 (IPesert et all 1201 lal) . 
To derive the quantities we have interpolated the BaSeL spec- 
tra at temperatures corresponding to the stellar effective temper- 
atures (r eff = 6000 K and T e $ = 5000 K for HD 209458 and 
HD 189733, respectively) and assuming AT = 1250 K for both 
cases. The Table also lists the measured planetary radius and the 
related uncertainty, as reported in the literature. 
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Fig. 5. Corrected Planetary Radii vs. Wavelength. Top panel, HD 
209458b. Middle panel, HD 189733b. Bottom Panel, GJ 436b. The open 
o circle indicates the planetary radius inferred from the literature, the 
x cross the corrected radius with f = 0.004 for HD 209458b and 
fo = 0.01 for HD 189733b and GJ 436b and the * asterisk the cor- 
rected radius with f = 0.03. The plotted points are shifted horizontally 
for an easier identification. 
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Fig. 6. Top panel: the synthetic light curve of the transit of C oRoT-2 as 
computed with the parameters in Table 1 of Alonso et al. ( 2008) (dashed 
line) together with the unperturbed profile obtained from Eq. ^ (dia- 
monds) and the corresponding best fit (solid line). The flux is measured 
in units of the out-of-transit flux level in all the cases. Bottom panel: the 
residuals of the best fit to the unperturbed transit profile in relative flux 
units vs. the orbital phase. 



In Fig. we show the corrected planetary radii vs. wave- 
length, along with the uncertainties reported in the literature . 
The a dopted stellar radii are ^hd 209458 = 1.14 6 R Q (Brow net all 
2001), Ri rojsgTO = 0.788 R Q dBaines et all l2008h . ™ = 
0A64R G dTorresl 120071) . respectively. The most active star (i.e. 
HD 189733) shows significant discrepancies between the planet 
radius inferred from the literature and the one determined here; 
this is an expected result since AR p oc / (see Eq.Q~2]in Sect. l2.21 > 
and may have interesting consequences if one wants to charac- 
terize the planet atmosphere through transmission spectroscopy. 
Similarly, in the quiet cases, even though our corrections are 
within the error bars reported in the literature (of the order of 
the hundreth or of the thousandth, see Table [5]), they may be- 
come significant on the determination of the physical properties 
of similar exoplanets. The wavelength dependence of the appar- 
ent radius of HD 189733b in the optical passbands has been ob- 
served by Sing et al.l d201 lbl) . They used an approach similar to 
that of Sect. 12.21 and concluded that the variation of Rp/R* vs. 
wavelength A cannot be reproduced by an equation analogue to 
Eq. dT3l ). thus providing evidence of a wavelength dependence 
of the planetary atmospheric absorption. 



3.3. Effects of stellar activity on orbital parameters 

To illustrate the impact of starspots on the determination of all 
the transit parameters, we cons ider the case of CoR oT-2 because 
it has been already studied by ICzesla et alj d2009l) and its spot 
activity has been investigated in detail. Since we are interested 
in the systematic errors produced by starspots, we simulate a 
noiseless transit profile adopting the t ime sa mpling and the pa- 
rameters as derived by Alon so et al.l d2008|). and apply E q. ([8]l 
with spot filling factors estimat ed from [ Lanza et all (120091) for 
the out-of-transit light curve and lSilva-Valio et a" for the 

transit light curve, respect ively. To derive the transit profile to 
be fitted with their model, lAlonso et al.1 (120081) averaged all the 
transits observed by CoRoT. Therefore, we assume that the spot 
perturbation corresponds to the average values of the filling fac- 
tors /o and f\ as derived from the modelling of the out-of-transit 
and the transit light curves, re spectively. By fittin g the light curve 
of CoRoT-2 outside transits. Lanza et a l. (2009) derived an av- 
erage spot filling factor = 0.07 adopting a spot contrast in 
the C oRoT white passband of Alanza = 0.665. On the other 
hand, ISilva-Valio et al.1 d2010h found an average filling factor 
of /i = 0.15 for a spot contrast of Asilva = 0.55. Multiplying 
it by the ra t io Al anzaMsilva we convert it into the scale of 
lLanza et a l. (2009), giving fi - 0.18. With those values of the 
filling factors and = 0.665, we compute the unperturbed tran- 
sit profile Do from Eq. ([8} and then fit it with the Levenberg- 
Marquardt algorithm as implemented in the IDL procedure lm- 
fit. The unperturbed transit profile and its best fit are plotted in 
Fig.|6]together with the initial profile assumed to be perturbed by 
the effects of the spots. The best fit has a reduced^ 2 = 1 .29 com- 
puted by assuming a relative standard deviation of the simulated 
da ta of 1.09 x 10~ 4 as in the case of the CoRoT data analysed 
bv I Alonso et al.l d2008l) . The best fit parameters with their stan- 
dard deviations as derived from the covarian ce matrix with the 
above assumed data standard deviation (see iPress et ail fl992, 
§15.4) are: a/R* = 6.76 ± 0.03, R p /R+ = 0.1726 ± 0.00006, 
i = 88°44 ± 0°18, m = 0.40 ± 0.04, and u 2 = 0.097 ± 0.04. 

We confirm that the radius variation is close to that estimated 
with the simple approach of Sect. 12.21 because the depth of the 
profile at mid transit gives the principal cons traint on the relativ e 
radius. However, contrary to the results of ICzesla et alJ (|2009), 
we find that the relative semimajor axis a/R* and the inclination 
i are also affected, with variations that exce ed three standard de - 
viations for both parameters (cf. Table 1 of lAlonso et all 2008). 
Note that the transit duration tj is the same in our model for both 
the unperturbed and the perturbed transit profiles. Therefore, the 
variations of a/R+ and i combine with each other as to main- 
tain tj practically constant (see Sect. l2.31 >. The variation of a/R* 
has an impact on the measurement of the density of th e star (cf. 
ISeager & Mallen-Ornelasi 120031: ISozzetti et all l2007h that, to- 
gether with its effective temperature, is used to derive its posi- 
tion on the H-R diagram and hence its age from the isochrones. 
We shall not consider this application any further because it is 
beyond the scope of the present work, but conclude that starspot 
effects cannot be neglected in the derivation of the stellar density 
from the transit fitting in the case of active planetary hosts. 

The limb-darkening parameters are not significantly af- 
fected. The small systematic deviations in the ingress and egress 
of the profile do not lead to any significant adjustment of the 
limb darkening (cf. upper and lower panels of Fig. |6j. 

Other stellar parameters can be affected by starspots, notably 
the effective temperature T e ff. The problem has been disc ussed 
in detail in the case of CoRoT-2 bv lGuillot & Havel d20TTb who 
found that the effect of spots is that of increasing the uncertainty 
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in the determination of T e g because of their changing covering 
factor along the activity cycle of the star. We shall not consider 
further these aspects because here we focus on the parameters 
directly related to the modelling of the transit light curve. 

4. Conclusions 

The observation of moderately active stars simultaneously in the 
optical and in the near infrared passbands allows us to discrim- 
inate between planetary transits and activity-induced variations, 
unless starspots are 2000 - 2500 K cooler than the unperturbed 
photosphere, which is generally not observed in main-sequence 
stars. This approach can be applied to individual light curve dips 
observed with a signal-to-noise ratio of* 10 - 15% of the flux 
variations in two well-separated passbands, e.g., U and K, and 
does not require a characterization of the stellar microvariability 
timescales which implies observations extended for at least 2-3 
stellar rotations. 

Moreover, infrared observations are important for the deriva- 
tion of the planetary radius and accordingly of the planetary den- 
sity. The correction of the effects due to starspots is necessary to 
avoid systematic errors in trasmission spectroscopy, especially 
when comparing the planet radius at optical and infrared pass- 
bands. A correct determination of the planetary radius allows 
us a better study of the planetary atmosphere and its composi- 
tions in the framework of the transmission and reflection spec- 
troscopy. We provide simple formulae, tabulations, and graphics 
to help estimating the starspot perturbation in the case of typical 
transiting systems. Moreover, we find that starspots may affect 
the relative semimajor axis a/R^ and the inclination i as derived 
by transit fitting. This may have a non-negligible impact on the 
determination of the stellar density, especially when optical light 
curves are used. Near infrared light curves are less affected, but 
the effect may still be significant for very active stars. 

A simultaneous broad-band photometric observation, from 
visible to infrared wavelengths, of a significant sample of stars, 
on one hand might reveal itself as a relevant tool to identifying 
and characterizing their stellar activity, allowing us to disenta- 
gle its contribution from planetary transit signals, on the other 
hand would help us to derive relevant and more accurate infor- 
mation on the primary pa rameters of the planets (as already sug- 
gested bv lJha et aLC OOO) and the physical properties of their at- 
mospheres, such as better constraints on temperatures and com- 
positions. 
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